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Abstract 



Non-abelian (NA) topological string defects exist in QCD with the spon- 
taneous chiral symmetry breaking U(Nf)L x U(Nf)n — * U(Nf)y. Anomaly 
effects lead to domain walls connected to these strings. In the framework of 
linear sigma model we find the configuration of these defects. Also we show 
that in this model metastable CP violating vacuua exist. Strings as well as 
extended regions of metastable vacuua may form during the chiral transition 
which could have interesting effects. 
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I. INTRODUCTION 



QCD with Nf flavors of massless quarks possesses the chiral symmetry U(Nf)i, x U(Nf)n. 
The generators of its diagonal subgroup U(Nf) v are parity invariant while the parity-odd 
generators give the axial transformations U(Nf)j L . (The latter do not form a group.) For low 
energy and temperatures the axial vector part of this symmetry suffers spontaneous symme- 
try breakdown (SSB), this is indicated by U{Nf)y x U(Nf)A — > U(Nf)v [0. This creates 
Nj massless pseudoscalar Goldstone bosons known as mesons. However in the real world, 
chiral symmetry is explicitly broken because of non-zero quark masses. This makes the 
pseudoscalar mesons massive. Even in the chiral limit, anomaly effects break the SU(Nf)y- 
invariant axial U(l) subgroup U(1)a C U(Nf)A to Z^ f p making the 77' meson heavier 
compared to other mesons. Finite temperature lattice QCD studies show that chiral sym- 
metry (but for U(1)a of course which is only effectively restored) is restored at temperatures, 
T larger than a critical value T x . Such chirally symmetric phase of matter is believed to 
have existed in the early Universe. The matter in the initial stages of the fire ball produced 
in heavy ion collisions is likely to be in the chirally symmetric phase. To see such a phase 
transition in laboratory is one of the motivations of heavy-ion collision experiments. 

The chirally symmetric and broken phases of QCD are characterized by an order param- 
eter (OP) which is the quark-antiquark condensate known as the chiral condensate $. It is 
an Nf x Nf matrix at each x. In the chiral symmetry limit (with no anomaly), all possible 
values of <3> constitute the order parameter space (OPS). In the broken symmetry phase in- 
stead, it is diffeomorphic to U(Nf). U(Nf) as the OPS allows for the existence of "abelian" 
as well as "non-abelian" topological string defects. Brandenberger et al. [f§ first studied the 
formation of "abelian" strings due to SSB of the U(1)a part of chiral symmetry. In || it 
was argued that these string defects will exist as long as the anomaly effects are small and 
will decay as the anomaly becomes substantial. In a previous work we have considered the 
effects of anomaly on abelian strings Q. We have argued that abelian strings will exist in 
the presence of anomaly with the typical structure of Nf domain walls connected to the 
string. Now, a closer look at the full OPS=U(Nf) shows that an additional class of topo- 
logical defects, known as non-abelian (NA) strings can also exist, as the arguments in ref. 
|| for example show. Note that here we use the phrase "non-abelian strings" in same way 
as the phrase "non-abelian monopoles" |5|U are used and not in the sense that ttx(OPS) is 
non-abelian. Here we find and discuss numerical solutions for these NA strings and note the 
important differences between the NA and abelian strings. Also we show that in the linear 
sigma model we employ, metastable vacua exist at zero temperature. Such metastable states 
have been found to exist in the non-linear sigma model |7j] before, but for temperatures close 
to the deconfinement and/or chiral phase transition. 

In the next section we briefly discuss how the non-abelian string defects arise from the 
topological structure of the OPS=U(Nf) during the SSB of chiral symmetry. In section III 
we find the numerical profiles for the NA strings using the linear sigma model Lagrangian 
M. In sec. IV we show existence of metastable CP violating states in the above model. 
In section V we conclude with some speculative remarks on the formation of these string 
defects and metastable states during and/or after the chiral phase transition. 
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II. NON-ABELIAN STRINGS IN QCD 



Topological defects are usually formed in phase transitions associated with SSB with 
non-trivial OPS. The type of defects formed depends on the dimension n of physical space 
and the topology of the OPS. When the OPS after SSB is discrete, the topological defects 
are domain walls in 3-dimensional (3-D) physical space. When the OPS after SSB is a circle 
the defects are vortices in 2-D and strings in 3-D physical spaces. Monopoles arise in 3-D 
when the OPS after SSB is a 2-sphere (S 2 ). For an introduction to and review of topological 
defects, see [[|. For our discussion we need to know only about the string defects. In the 
following we discuss briefly about how string defects arise. 

Topological string defects arise when there exist non-trivial loops in the OPS, for example 
when it is a circle. A loop going around the circle k times cannot be smoothly shrunk to a 
point or cannot be smoothly deformed to another loop going I ^ k times around the circle. 
So there exist non-trivial loops, characterized by windings, on the circle. There is a close 
connection between the non-trivial loops and the string defects in physical space, see Fig.l. 
A string defect of winding k corresponds to a non-trivial loop winding k times around a 
circle, not contractible to a point in the OPS, as a loop 1 around the string is traversed once 
in physical space. In Fig. 1 we show an example of winding number 1 string defect for the 
case of a complex scalar field \l/ with real and imaginary parts and ^2 when SSB is of 
U(l). The figure on the left of Fig.l shows the typical effective potential for such a case. 
The minima M of the potential is a circle with radius The configuration \I/ of the string 
on the loop 1 (on xy-plane) around the string core is shown on the right. Magnitude of \I> 
is represented by the length of the vectors and the phase by the angle 9 these vectors make 
with the positive x-axis. As one goes around 1 once \l/ on it traces a complete winding on M. 
However, the configuration of \I/ could have been such that it winds fc-times (for a winding 
number k string) around M for a single loop on 1. 

The stability of the string defect is purely due to topological reasons. The total variation 
A9 of 9 around 1 can have only discrete values, namely integral multiples 2irk of 2tt. So by 
smooth changes of \1/ on 1, one can not change A9. Also by smoothly shrinking or enlarging 
1, such that it does not cross \I> = when 9 becomes undefined, one cannot change A9 if 
the underlying \1/ is continuous. On the other hand if 1 can be shrunk to a point ( keeping 
\l/ continuous and non-zero except perhaps at that point, and A9 fixed in the process), then 
either k — or \& must vanish at that point as it is continuous. By definition, k 7^ when 
there is a string (defect). The point where \I/ = when k 7^ is the core of the string. 
Strings in 3-D are either infinitely long or exist in loops. As strings have non-zero string 
tension, a string loop is unstable to shrinking. 

It is important to note that string defects arise because one can have loops with non- 
trivial winding on M which is a circle in the above example. Such non-trivial loops exist in 
U(Nf) which is the OPS after SSB of chiral symmetry. Thus let T a (a e{l, Nj — 1}) be the 
SU(Nf) Lie algebra generators (the analogues of the Gell-mann matrices A a ) and 
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FIG. 1. On the left is the effective potential for SSB of U(l). The vacua are represented by 
the dashed circle. On the right is shown the configuration of the complex scalar field on the loop 
1 around the string core. 
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a particular non-abelian string configuration is 
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where 9 is the angle on 1. The first factor is an open curve in SU(Nf) from 1 to its center 

i— -i — 

e N f 1. The second is an open curve in the center £7(1) of U(Nf) from 1 to e N f 1. On taking 
their product as in the above equation, the values at 9 = 2ir multiply to 1 giving a closed 
string. The powers M(9) N (N 6 Z\{0}) and their (continuous) deformations produce all 
possible strings, the deformations not changing their topological types. The beginning and 
end points of the non-abelian and abelian factors remain fixed during these deformations. 
Note that in M(9) nN f (n G Z\{0}), the non-abelian factor 



inO 



N f (N f -l) 
2 



(4) 



is a closed loop in SU(Nf). Such loops can be deformed to a point as the first homotopy 
group of SU(Nf), iri(SU(Nf)) is 0. Thus M(9) nNf gives the abelian strings. In other words, 
the elementary strings are given by M(9) and M(9)- 1 . Note that by construction M(9) n 
gives a loop winding n times. Consider the loop given by 
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for which the open curve in SU(Nf) remains the same as that of M(6), but one traverses 
the two end points in U(l) in an anti-clockwise direction. This loop is actually homotopic 
to the loop given by M(8) 1 ~ Nf . That is to say, there is a string loop homotopic to that of 
M(9) 1 ~ N f which gives the string configuration of (4). 

In the above we have considered SSB of chiral symmetry However in the real world 
this symmetry is explicitly broken due to quark masses, and also the U(1)a part of U(Nf)A 
is broken down to Z^ f due to axial anomaly and instanton effects. The effect of anomaly 
changes the usual cylindrical structure of the string creating domain walls connecting to 
the string core. In the case of abelian strings, for Nf = 3, the string gets connected to 
three domain walls [f|. In the case of elementary NA strings the number of domain walls 
can be either one or two as we explain in the next section. The effect of quark masses 
makes the string configurations non-static. However the string defects are stable against 
small fluctuations and can only decay due to annihilation with other defects with opposite 
windings. 

In the next section we will numerically find the configuration of the elementary non- 
abelian strings for the case of Nf = 3. For our calculations we use the linear sigma model 
Lagrangian considered in ref. ||. Numerically it is easier to find the static configurations. 
Because of this we work in the chiral limit which allows for static solutions. 



III. CONFIGURATIONS OF THE NON- ABELIAN STRING 

To find the field configurations for the non-abelian strings we consider the following 
U(3)l x U(3)r linear sigma model for 3 quark flavors fiJlCD, 



£($) = Tr (d^d^ - m 2 $ f $) - Ai [Tr (V$)] 2 - A 2 Tr (V$ 



+ c 



Det ($) + Det (V)] + Tr [#($ + $ f ) 



(6) 



where we take c > for specificity. $ is a complex 3x3 matrix of fields of the scalar and 
pseudoscalar mesons: 



$ = T a (p a = T a (a a + in a ). 
The field $ transforms under the chiral transformation U(3)l X U(3)r as 
Q^UtQUl , Ui, r = exp (i w% r T a 



(7) 



(8) 



One can rewrite the infinitesimal left-right transformations in terms of vector-axial vector 
transformations with parameters Wy A = (wf ±u>")/2. The infinitesimal form of the U(3)l X 
U(3)r symmetry transformation (H) is 



$_^$ + ^ [T a ^]+iu a A {T a ,$} . 



(9) 



$ is a singlet under the U(l)y transformations associated with exp(iw v T°). In QCD, they 
give rise to a conserved charge identified with the baryon number. 

In the above model the determinant term takes into account the instanton effect which 
explicitly breaks the U(1)a symmetry [|K|]. It is not clear whether such a term can be justified 



for high temperatures, because vanishing <£> makes the instanton effects to disappear ||11|| , 
though lattice results show that U(1)a is effectively restored at slightly higher temperatures 
IJT2]] than the temperature at which SU(N)a symmetry is restored. The last term where 
(H = T a h a , h a = constants) is due to non-zero quark masses. When c — H — and Ai > 0, 
and A2 > 0, the Lagrangian has a global symmetry U{Nf)y x U(Nf)A for m 2 > 0. For 
m? < the axial symmetry U(Nf)A is spontaneously broken to identity and the OPS is 
U(Nf). This results in Nj Goldstone bosons. For Nf = 3 these Goldstone bosons are the 
7r's, K's, t] and rf . However when just c 7^ 0, the U(1)a is further broken to Z% by the axial 
anomaly. U(3)a symmetry is, in addition, explicitly broken by non-zero quark masses. Also 
SU(3)y is explicitly broken due to the difference in u, d, and s quark masses giving rise to 
splitting in the meson masses. 

In order to find the static non-abelian string configurations we take H to be zero. Even 
when H = 0, not all the non-abelian strings are static in the presence of anomaly just like 
the abelian string. To illustrate this we consider the NA string corresponding to the non- 
trivial loop M(9) [with N f = 3] on U(3). Scaling M(9) by the vacuum $ = ^1, we get 
the loop in the space of fields $ to be 

1 ( a ° 

$(0) = -±M{6) = — (T I . (10) 



x/6 VQ 



(Toe- 



As only the last diagonal component of $ varies with 9, the anomaly contribution to the 
effective Lagrangian from such field variation is 

alcosO. (11) 



3^6 

For such a term, the ground state prefers the unique value cosO = 1. Because of this, the 
9 distribution around the string will not be symmetric. In most of the region around the 
string, 9 will take the value which minimizes the potential energy term (10). Non-zero 9 
values will be confined to a wall-like region. Because of non-zero surface tension of the wall, 
motion of the string towards the wall decreases the energy of the configuration. So one will 
not have static solutions in this case. 

On the other hand, the loop corresponding to M(9) n in OPS gives rise to a string 
configuration with |n| domain walls. This configuration is static as the motion of the string 
along any domain wall does not decrease the energy of the configuration. To find out the 
details of one such NA string configuration, we consider Nf = 3 and n = —2. Here M(9)~ 2 
gives a loop homotopic to the loop of M'(9). For Nf = 3, M'(9) is given by exp(i9T s ) where 
T a is 



T s = 1 . (12) 
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The corresponding loop in the space of fields $ is from 




m = ^=M'(9) = -L I o roe* . (13) 
Writing \/§(0i + ifo) for <7oe td , the effective Lagrangian for such a $ deduced from (6), 



with Ai = A, A2 = and H = 0, is given by 



where 



C = ^01^01 + ^02^02 - V(0!, 2 ), (14) 



V(01, 02) = (0? + $ + \ (0? + 0^) 2 

The values of the parameters in the above equation are A = 53.5322, c = 6685.19MCV and 
m 2 = — (591.256MeV) 2 . For these values of the parameters we get a = 110.23MeK, the 
sigma meson mass m CT = lOOOMeK, and the 7/ mass m v i = 950MeV. Note that the anomaly 
term in V expressed in 9 is — ca /2\^Qcos(29) for which 9 = and 9 = n are degenerate. 
This leads to a string configuration attached to 2 domain walls. In the following we will find 
the approximate numerical profile of the string associated with M'(9). 

Considering the string along the ^-direction, 0j for the static string satisfy the following 
field equations: 

2 dyj^M 
v 01 " ' 

V 2 02 = dV f^\ (16) 



where V 2 = + Jp. Since in the presence of anomaly, the string profile is non-cylindrical, 
one has to find the solution by an energy minimization technique. For this we consider an 
initial cylindrical configuration with 0« = at t = and evolve it with the field equations 
derived from Eq.(6) with a dissipative term: 

d 2 h + a #i _ = _5V(0i 1 02) 

rft 2 dt * ! (90i 

Gp02 rf0 2 2 9V(0l,0 2 ) , , 

— + a — -V0 2 = ^— . (17) 

Here a is the dissipation coefficient. The dissipation term has been used to converge the 
initial configuration (the circular configuration, which is not the correct configuration in 
the presence of the anomaly) to the approximate configuration of the string connected to 
domain walls. Once the evolution is done for a sufficiently long time, the configuration 
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FIG. 2. Configuration of the string for c = 6685.19Mey and H = 0. Figure on the left gives 
— \J4>1 + $2 = —<p an d figure on the right gives the vector field plot of the full field 0i(r) + ifair)- 

does not evolve with time even after the dissipation is switched off, which suggests that 
the configuration is more or less static and so an approximate solution. On evolving the 
initial cylindrically symmetric configuration with the above equations, the domain walls 
connected to the string develop. In the numerical work, we considered the n' mass to be 
the experimental value m r f ~ 950MeV and set other pseudoscalar meson masses to zero. 
In FIG. 2 we show the configuration of the string with the domain walls joined to it. In the 

left figure we plot —\J4>\ + 02, and in the right plot the vector field of $(r) = (0i(r), ^(r)). 
The magnitude of the field is proportional to the length of the vectors and the angle these 
vectors make with +ve x-axis is the phase. Clearly there is no symmetric distribution of 
the phase around the string. In the core of the string the full chiral symmetry associated 
with T s is restored. That is because in the string core, as we explain below, $ becomes 
diag(0, 0, constant) which is invariant under transformations generated by T s . Although the 
field configuration has the structure of domain wall, the symmetry is restored only in the 
core of the string. It is not clear if it is the string tension or the surface tension of the domain 
wall that will dominate the dynamics when these configurations are formed after the chiral 
transition. 

In the absence of anomaly the NA and abelian string configurations j|] look superficially 
similar. But there are crucial differences. For in the core of the abelian string all the 
components of $ vanish and so the full chiral symmetry is restored. This is because the 
phase of all components of $ change by 27m (n being a non-zero integer) as one goes round 
an abelian string obliging all these components to vanish at the core of the string. But 
for a NA string one or more components of $ need not change in phase as the string is 
encircled, and hence need not vanish in the string core. As a consequence the full chiral 
symmetry associated with a particular linear combination of T and T N 2_ 1 for which the 
nonzero components are singlets is restored. The other difference is as we have mentioned 
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above, for c ^ 0, the elementary abelian string for Nf flavours is connected to Nf domain 
walls while the elementary NA strings (associated with M(9) and M(9)~ 1 ) are connected to 
1, 2, Nf — 1 domain walls. (But note that M'(6) has two domain walls.) 



IV. METASTABLE STATES IN 17(3) x f7(3) LINEAR SIGMA MODEL 

As we have mentioned, in the chiral limit, anomaly breaks U (1)^ to Z3 in 3 flavour QCD. 
This leads to a set of discrete ground states related by discrete transformations of C7 (3). It is 
simple to find out the ground states in the chiral limit. Let us consider the effective potential 
in Eq.5, 

V($) = Tr (m 2 + Ai for 2 + A 2 Tr ($t$) 2 

-c[Det($) + Det(V)] - Tr [i7(<l> + $ f )] . (18) 

In the chiral limit when H = 0, one can choose the ground states to be of the form $ — 
oqTq/\/2. When the anomaly is also absent, c = 0, and all U(3)a rotations of $ cos t 
no energy. For c ^ 0, only a few discrete axial rotations do not change the potential 
energy. Thus there are a few ground states. Some of the ground states are obtained by the 
discrete rotations U&oU, where U = exp(—i9i/2), with 9\ = ^r, 4^. Also one can choose 
U = exp{—i9 2 T /2), where T is a linear combinations of T a 's and T , for example T s in Eq.ll. 
For T = T s , 62 = 7t~ gives another ground state. All these ground states, except $0, will have 
non-zero U{1)a phase giving rise to CP violating effects f7j. 

For the full potential (18), these states with 9x t2 7^ will be of higher energy. Still 
these states with differing values of 9 12 are separated by energy barriers. When the quark 
mass is switched on, only one of these states becomes the true ground state. However, 
for small quark masses, the barrier between the ground state and other states may remain 
appreciable making them metastable. Because of their CP violating nature it is interesting 
to check if such metastable states can exist after including the effects of the meson masses 
and their splittings. In the non-linear sigma model such states are known to appear 
for temperatures close to the deconfinement transition. Here we now show that such states 
exist also in the linear sigma model at zero temperature. Since these are due to effects of 
anomaly we expect that they will survive as long as anomaly effects are significant. 

In order to find metastable states, we took the values of the parameters in Eq. (18) 
from || which give an approximate description of meson masses and their splittings. The 
ground state is <3> = a T /\/2 + cfsTg, where a$ ~ ISOMeV and a 8 ~ —2AMeV while H is 
h T /V2 + h 8 T 8 , where h = (286.094Me\/) 3 and h 8 = -(UttSmMeVf . Finally \ x = 1.4 
and A2 = 46.484. In order to see possible metastable states, we fix a and cr 8 and plot the 
#1,2 dependence of the potential energy. In case there is a local minimum for non-zero 9^2, 
we then minimize V($) relaxing <f and a 8 fixing # 1)2 . If the new <jq and a 8 are close to their 
initial values, then it implies the existence of metastable states. We find no metastable states 
for non-zero 9 X , so there exists no metastable state with purely U(1)a phase. For # 2 there 
exists a metastable state at 9 2 = tt. In Fig. 3 we plot the 9 2 dependence of the potential 
energy for m = and m q 7^ 0. The figure clearly shows that there is a local minimum at 
9 2 = tt for m q ^ 0. Fixing # 2 at ir we minimize V($) varying a Q and <f 8 . The new values are 
cfo ~ 120MeV and <f 8 = —28MeV indeed very close to the initial values. 
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FIG. 3. #2 dependence of the potential energy for m q = and m q ^ 0. 

CONCLUSIONS 

In summary, we have numerically found an approximate static solution of a non-abelian 
string (domain wall) in the SSB of U(Nf) x U (Nf) — ► U(Nf) in the linear sigma model. Also 
we have shown that metastable states exist in this model for realistic quark masses. Abelian 
and non-abelian extended regions in metastable vacua may form during the chiral transition 
ocurring in the early stages of heavy-ion collisions as well as in the early Universe. However, 
formation of strings and large metastable domains are complementary as formation of the 
one suppresses the formation of other. More the number of uncorrelated domains, more are 
the number of strings, with less probability for an extended domain in a metastable state. 
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